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The QCD phase diagram at finite temperature and density is a topic of considerable interest. 
Ahhough much progress has been made in recent years, some open questions remain. Even 
at zero density, the order of the transition for two light flavors of fermions has not yet been 
conclusively established. While considerable evidence exists in favor of a second-order transition 
for massless quarks and a crossover for massive quarks, some recent results with two flavors of 
staggered fermions suggest a transition of first order 

Since lattice simulations are performed in finite simulation volumes, actual phase transitions can- 
not be observed directly. Thus, finite-size scaling is a very useful tool in the analysis of lattice 
data. By comparing the scaling behavior of observables to the expected scaling properties, values 
of critical exponents can be confirmed and the order as well as the universality class of a transi- 
tion can be established. In the comparison to lattice QCD results, the critical exponents and the 
universal scaling functions have been obtained mainly by means of lattice simulations of 0(N) 
spin models, and results are usually restricted to the critical temperature or the point at which the 
susceptibilities peak. 

We propose to use a non-perturbative RenormaUzation Group method for this purpose. We have 
calculated the critical finite-size scaling behavior and the universal scaling functions for the three- 
dimensional 0(4)-model for a wide range of temperatures and values of the symmetry breaking 
parameter. Our results are suitable for a comparison to lattice QCD results for the chiral suscep- 
tibility and the order parameter and can be used to check the consistency of the finite-size scaling 
behavior with that of the 0(N) universaUty class. 
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1. Introduction 

In QCD, two phase transitions take place at finite temperature and density, a deconfinement 
phase transition dominated by the gauge fields, and a chiral symmetry restoration transition. Whether 
these transitions coincide is still not settled [|l|, If the QCD phase transition is dominated by 
the restoration of chiral symmetry, the conventional scenario [Q] holds that the transition is of sec- 
ond order for Nf = 2 massless quark flavors, and a crossover for massive quarks. While there is 
much evidence in favor of this scenario [^], results obtained with two dynamical flavors of stag- 
gered fermions suggest a first-order transition |^, and recent work by the same group supports 
these findings |^, ^. 

Since lattice simulations are necessarily performed in finite simulation volumes, and since 
strictly speaking phase transitions only occur in the infinite-volume limit, it is difficult to conclu- 
sively determine the order of a phase transition. A further complication is the explicit breaking of 
the chiral symmetry by the finite quark mass, which corresponds to an external symmetry-breaking 
field. An important tool is therefore the finite-size scaling analysis of the simulation data. Due 
to universality, in the case of critical behavior results are expected to fall onto universal scaling 
curves, characterized by critical exponents and scaling functions of a particular universality class. 

Necessary for establishing the nature of critical behavior in a finite volume is therefore knowl- 
edge of two things: of the critical exponents, and of the scaling behavior. So far, scaling functions 
have mainly been determined from lattice simulations of 0(N) spin models |^[T^]. While there is 
evidence of scaling for Wilson fermions, the expected infinite-volume scaling behavior is not seen 



for staggered fermions [11, 12]. Results with a modified QCD action with two flavors of staggered 
fermions suggest that current simulation volumes might be outside the scaling region [pj|]. 

In the work presented, we use a non-perturbative Renormalization Group (RG) method to 
calculate critical exponents as well as the finite-size scaling functions. While an abundance of 
results exists on critical exponents, this is to our knowledge the first time that finite-size scaling 
functions have been obtained with non-perturbative RG methods. As an advantage, a very wide 
range of parameters and volume sizes is accessible to the method, and thus direct comparisons to 
lattice simulations are possible and the size of the scaling region can be assessed. 



2. Renormalization Group method 

We use a non-perturbative Renormalization Group method to calculate observables such as 
the order parameter and the susceptibility. The method accounts for long-range fluctuations and is 
suitable to describe critical behavior. For a review of functional RG methods, see e.g. [14, 15, p^ ]. 



In the framework of the effective action, an infrared cutoff k is introduced, and by varying the 
cutoff scale quantum fluctuations are integrated out in a systematic way. The values of coupling 
constants at an initial UV cutoff scale A are the input. After the IR cutoff has been lowered to 
k = Q, all quantum fluctuations below A have been integrated out and are included in the values of 
the couplings in the effective action. The change of the couplings under a variation of the cutoff 
scale k is governed by a flow equation. For a particular choice of the RG cutoff scheme, the flow 
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p V T] d 

J. Engels et al. p] lattice 0.380 0.7423 (0.024) 4.86 



D.F. LititmandJ. M. Pawlowski[18] RG (0.4022) 0.8043 - (5.00) 



this work RG 0.4030(3) 0.8053(6) (0.0046(4)) 4.9727(5) 

Table 1: Critical exponents for 0(4) in c/ = 3. Values in brackets are calculated using scaling laws. 



equation for the i/-dimensional 0(N)-model in infinite volume reads [17, 18, 19] 



^^[/,(a,7r)- ^^^^^^2+1) \k^+Ml{k)^ k^+Ml{k) J ' ^^"^^ 

We have adapted this cutoff scheme to calculations in a finite volume [pO|], and are currently inves- 
tigating the dependence of our finite-volume results on the choice of cutoff scheme. The effective 
potential is expanded in local «-point couplings around its minimum (Jo{k) 

Uk{o, 71) = ao{k) + ai {k){o'^ + tt^ _ aQ{kf) + a2{k){o'^ + tt^ _ OQ{kff + ...-Ha, 

where H is the fixed, external symmetry -breaking field, and the additional condition 2a i (k)ao{k) = 
H keeps the minimum at (a, n) = (ao(^),0). The RG flow equation is solved numerically. 

In <i = 3 dimensions, the initial value of the minimum at the cutoff scale ao(A) serves as proxy 
for the temperature, (do (A) — aQ"°'^^'(A)) ~ (r — Tc). For this first calculation, we have chosen 
a cutoff scale A = 1.0 GeV, which is comparable to a typical lattice cutoff {7i/a « 1.5 GeV) in 
current lattice simulations. 



3. Scaling in Infinite Volume 

Close to a critical point, such as a second-order phase transition, the correlation length ^ di- 
verges, and long-range fluctuations dominate the behavior of a system. As a consequence, certain 
observables become independent of the short-scale details of the system. Close to the transition, 
their behavior is described by only a few critical exponents and scaling functions which are char- 
acteristic of the universality class of the transition. For example, the order parameter behaves 
according to 

M{t,h)=h'^^f{z), z = t/h'l\ (3.1) 

where f{z) is the scaling function and z is the scaling variable. The two conditions M{t,h = 0) = 
{—t)^ and /(O) = 1 determine normalization constants for the dimensionless temperature and field 
parameters, t = {T — Tc)/To and h = H/Hq. 

We obtain values for the critical exponents j8, V, and 5 by fitting directly to the observables 
M = 00 and ^ = I /M^ as functions of the temperature and the symmetry-breaking field. In Tab. |l], 
our results are compared to those of a different RG calculation, where we find perfect agreement, 
and lattice spin model calculations, where the agreement is reasonably good and the deviations can 
Ukely be explained by the restriction to local couplings in our approach. 

In Figs. |T] and |^ we validate the results for the critical exponents by confirming the scaling 
behavior. In Fig. [T],the order parameter M is shown as a function of the temperature t for several 



3 



Finite Size Scaling for the 0(N) universality class from Renormalization Group Methods Bertram Klein 




M/K 



1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 

H /MeW^'^ 






— 100.0 






— 90.0 






— 80.0 






— 70.0 






— 60.0 






— 50.0 






— 40.0 






— 30.0 






— 20.0 






— 10.0 




: 0(4), f/ = 3 

■ .... 1 .... 1 ... . 


1 . ^-^-^-t— . . . . 1 





-30 



-20 



-10 







10 



20 



30 



Figure 1: Scaling for the order parameter M in infinite volume. M is shown as a function of the temperature 
t for several different values of the symmetry breaking field H (left panel). After rescaling M/h^l^ as a 
function of z, results for these values of H collapse onto a single curve (right panel). 
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Figure 2: Scaling for the susceptibility x in infinite volume, x is shown as a function of the temperature f 
for several different values of the field H (left panel). After rescaling x/h^'^^^^^^ as a function of z, results 
for these values of H collapse onto a single curve (right panel). 

different values of the field H (right). Plotting M/h^l^ against the scaling variable z, the different 
curves collapse onto the universal scaling function f{z) (left). In the second figure, scaling is 
validated for the susceptibility, again the curves collapse after rescaling. We have checked that 
the scaling behavior holds over three orders of magnitude for the field, from // = 1.0 to 1.0 x 10^ 
MeV^/^. This confirms that our approach incorporates the critical fluctuations responsible for the 
universal behavior. 

4. Finite- Size Scaling 

Since a finite volume provides a natural IR cutoff, the long-range fluctuations are affected and 
the critical behavior is changed when a system is put into a finite volume. Because these cutoff 
effects depend on the relative size of the correlation length t, and the volume size L, it is natural 
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Figure 3: Finite-size scaling for the order parameter M at z = 0. M is shown as a function of h for different 
volume sizes (left). The finite-size scaled order parameter ML^I^ is plotted against the finite-size scaling 
variable hd^^l" (right). 



to assume that deviations from the infinite-volume scaling behavior are controlled by this ratio. 



According to the finite-size scaling hypothesis [21], the ratio of a thermodynamic observable in the 



finite-size and the infinite-size system at the same temperature is a function of only this ratio t,/L. 
For example, according to this hypothesis one expects for the order parameter in finite and infinite 
volume 

"'<"-^fi«V (4.,) 



M^{t) 

If this hypothesis holds, universal finite-size scaling functions can be obtained by realizing that one 
needs to change t ~ L^^l^ with the volume to keep the ratio t, /L constant, since E, ~ To keep 
the scaling variable z. = t/h^^^^^^ constant as well, one also needs to vary h ~ L-l^^^". Taking the 
infinite-volume scaling behavior M{t,h) = h^^^f{z) into account, one expects the combination 

MLP^/^ = QM{z,hL^^^'') (4.2) 

to be a universal function for a given value of z. This is indeed borne out by our results. In Fig. ||, 
the deviation from the infinite- volume behavior for the order parameter as a function of h is readily 
apparent (left). For large values of h, where the correlation length is small, the results for different 
volume sizes all approach the infinite-volume limit. After the finite-size rescaling, the curves for 
the larger volume sizes once again are on top of each other (right), which clearly demonstrates the 
scaling behavior. For the susceptibility x, the finite-size scahng behavior is shown in Fig. For 
volume sizes 10 fm < L < 100 fm, the unsealed susceptibilities differ by two orders of magnitude, 
but agree within 10% after rescaling. 

For small volume sizes (for our choice of parameters L < 10 fm), deviations from the leading- 
order scaling behavior become apparent. These deviations can be explained from an RG analysis 
of finite-size scaling beyond leading order, which takes irrelevant operators into account. Eq. ( p!2[ ) 
then reads more fully with these corrections 

MlP/^ = dSiz^hL^'"") + ^Qi\z,hLf''^n + . • • , (4.3) 
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Figure 4: Finite-size scaling for the susceptibility x z = 0- X plotted vs. h for different volume sizes 
(left). The finite-size scaled susceptibiHty L'^^^X is plotted vs. the finite-size scaUng variable hL^^I^ (right). 



where the exponent CO is associated with the first irrelevant operator at the critical point. From a 
global fit to our results, we simultaneously determine the coefficient functions and the exponent ft). 
Our result (O = 0.74(4) is in good agreement with (O = 0.7338, which was obtained in [18] from 
an RG fixed-point analysis carried out in infinite volume in the same RG scheme. The results for 
the coefficient functions {z,hL^^/^) and {z,hL^^I^) determined from the volume range 10 
fm < L < 100 fm are quite robust and almost unchanged when the analysis is extended to also take 
the region down to L = 1 fm into account, even though for our current choice of parameters the 
deviations from the leading-order scaling are already quite sizeable below 10 fm. 

An important conclusion from this is that for small volumes, sub-leading corrections to the 
scaling behavior can become so large that the leading-order scaling analysis is no longer sufficient. 
While these results have presently been obtained only for a particular choice of parameters and 
the volume sizes cannot be immediately equated with those of QCD lattice calculations, they still 
suggest that scaling deviations play a role for small volumes in current lattice simulations. 



5. Conclusions 

In the work presented, we have investigated finite-size scaling in the 0(N) universality class 
with non-perturbative Renormalization Group methods. For the 0(4) model in <i = 3, the critical 
exponents and the scaling functions have been obtained, and the infinite-volume scaling behavior 
has been used to validate these results. 

Adopting the RG scheme to calculations in finite volume, we have further demonstrated that 
finite-size scaling behavior can be obtained in this framework. We have calculated finite-size scal- 
ing functions for the order parameter and the susceptibility for a very wide range of values for the 
temperature and the symmetry-breaking field and for a wide volume range. The deviations from the 
leading-order scaling behavior that we observe can be adequately explained by the scaling correc- 
tions expected from an RG analysis. However, since these corrections are sizable for small volume 
sizes, they counsel caution for an analysis that takes only the leading-order scaling behavior into 
account for lattice results from small volumes. 
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The results that we have obtained can be compared directly to results from QCD lattice sim- 
ulations and can be used to check compatibility with 0(4) scaling behavior. An extension to the 
0(2) universality class and such a direct comparison are in preparation. 
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